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Abstract. The fractional-order derivative of a non-constant periodic function is not pe- 
riodic with the same period. Consequently, any time-invariant fractional-order systems 
do not have a non-constant periodic solution. This property limits the applicability of 
fractional derivatives and makes it unfavorable to model periodic real phenomena. This 
article introduces a modification to the Caputo and Rieman-Liouville fractional-order 
operators by fixing their memory length and varying the lower terminal. It is shown 
that this modified definition of fractional derivative preserves the periodicity. There- 
fore, periodic solutions can be expected in fractional-order systems in terms of the new 
fractional derivative operator. To confirm this assertion, one investigates two examples, 
one linear system for which one gives an exact periodic solution by its analytical ex- 
pression and another nonlinear system for which one provides exact periodic solutions 
using qualitative and numerical methods. 
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1 Introduction 


The history of fractional calculus goes back to the end of the 17th century when L’Hopital 
asked Leibniz what meaning could be ascribed to D"f if n were a fraction? Since that, 
time-fractional calculus has drawn the attention of many famous mathematicians, such as 
Euler, Laplace, Fourier, Abel, Liouville, Riemann, and Laurent [18]. The advantages of 
fractional calculus have been described and pointed out in the last few decades by many 
authors [8, 15-19]. It has been shown that the fractional-order models of realistic systems 
are regularly more adequate than usually used integer-order models. Applications of these 
fractional-order models spread in many fields, such as viscoelastic systems, dielectric polar- 
ization, electrode-electrolyte polarization, electromagnetic waves, quantum evolution of com- 
plex systems, and so on [6,10,11,14,20]. There are three definitions most frequently used 
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for the general fractional differential operators. The first one is the Griinwald-Letnikov (GL) 
fractional differential operator defined by the limit of a fractional-order backward difference 
and has an advantage for numerical simulations. The second type is the Riemann-Liouville 
(RL) definition; this operator played a pivotal role in developing the fractional calculus the- 
ory. Using these two fractional differential operators in modeling real phenomena leads to 
mathematical models with initial conditions expressed in terms of fractional derivatives that 
do not have known physical interpretation. The third type is the Caputo derivative having the 
advantage of dealing models with initial conditions expressed in terms of the field variables 
and their integer-order derivatives, having clear physical interpretations [9]. Recently it has 
been demonstrated that the fractional-order derivative of a non-constant periodic function is 
not a periodic function with the same period [13, 22,23] and in [5] the authors studied quasi- 
periodic properties of fractional order integrals and derivatives of periodic functions. As a 
consequence of the non-periodicity of the fractional derivative of a T—periodic function, the 
time-invariant fractional-order systems do not have any non-constant exact periodic solution 
unless the lower terminal of the derivative is -oo [12, 13,23], which is not realistic. This prop- 
erty limits the applicability of the fractional derivative and makes it unfavorable for periodic 
real phenomena. In [1], the authors have proposed a modification of the Griinwald-Letnikov 
fractional differential operator, which consists of fixing the memory length and varying the 
lower terminal of the derivative. They have demonstrated that the modified definition of 
fractional derivative preserves the periodicity. The present paper extends this modification to 
the Caputo and Rieman-Liouville fractional-order operators. Tow examples are investigated 
to confirm that periodic solutions arise in fractional-order systems when the new fractional 
derivative operator is used. One linear system for which one gives an exact periodic solution 
defined by its analytical expression and another nonlinear system for which one provides an 
exact periodic solution using both qualitative and numerical methods. 


2 Fractional-Order Derivatives 


As said above, the most usual definitions of fractional-order derivative are the Griinwald- 
Letnikov, the Riemann-Liouville and the Caputo definitions [17]. For 0 < a ¢ IN, the a—th 


order derivative of a function f(t) with respect to f and a terminal value a is given in the sense 
of 


¢ Griinwald-Letnikov by 


DEF) = tim me (ay ( { ) fk), (21) 
=0 
nh=x-—a 


a \ — T(w+1) 
where ( k ) - ATEN" 


¢ Riemann-Liouville by 





DEO = reg maya C—O FaD). (2.2) 
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¢ Caputo by 
t 


5 / (¢— t)"—*1 f0" (yd, (2.3) 


a 


1 
m—k 





SDEF() = 5 


In (2.2) and (2.3), m is the first integer greater than a, and I’(.) is the Gamma function. 
The following theorems reveal a remarkable property for the fractional derivatives based on 
Caputo definition, Griinwald-Letnikov definition, Riemann-Liouville definition [22]. 


Theorem 2.1. Suppose that f(t) is a non constant periodic function with period T. 
If f(t) is m-times differentiable, then the functions ©D*f(t), where 0 < a ¢ IN and m is the first 
integer greater than x, cannot be a periodic functions with period T. 


Theorem 2.2. Suppose that f(t) is (m-1)-times continuously differentiable and f\")(t) is bounded. If 
f(t) is a non-constant periodic function with period T, then the functions C/D* f(t) and ®“D* f(t), 


a 
where 0 < « ¢ N and m is the first integer greater than x, cannot be periodic functions with period T. 


Example 2.3. Let f(t) = sin(t). One has 


Hence 


where 0 < a < 1 and Ey g(t) is the generalized Mittag-Leffler function defined by 
oo tk 
Fuel) = De Take BY 


k=0 


Numerical simulations showed that f*- "E9525 ( 17) is not a periodic function where 0 < a < 
1, even if a = 1 this function is the periodic function cos(f). 


As a consequence of the above theorems, periodic solution cannot be expected in fractional- 
order systems, under any circumstances [22,23]. 


Corollary 2.4. A differential equation of fractional-order in the form 
aDex(t) = F(x(t)), 
where 0 < « ¢ IN, cannot have any non-constant smooth periodic solution. 
This property highlights one of the basic differences between fractional-order derivative 
and integer-order one, and it makes fractional-order systems unfavourable for a wide range 


of real periodic phenomena. Therefore in this paper one overcomes this problem by imposing 
a simple modification to both Riemann-Liouville and Caputo definitions. 
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3. The Fractional-Order Derivative with Sliding Fixed Memory Length 


one first recalls the Griinwald-Letnikov fractional-order derivative with fixed memory length 
introduced in [1]. 


Definition 3.1. (The Griinwald-Letnikov fractional derivative with fixed memory length) 
Let a > 0, L > 0, m an integer such that m — 1 < a < m and f an integrable function in the 
interval [a — L,b]. The operator is D? defined by : 
L 
: T(a+1) 


“DESO = lin re OM are eae kh), t € [a,b], (3.1) 





is called the Griinwald-Letnikov fractional derivative with sliding fixed memory length. 


The following proposition gives an evaluation of the limit in the definition of Gritnwald- 
Letnikov fractional derivative with sliding fixed memory length. 


Proposition 3.2. Under the assumptions of definition (3.1), if the function f is m-differentiable with 
f™ €L[a—L,b], then 





m—1 ¢(k)(¢ _ k-a ; 
Mo De E(t) _ bs a = | os 5 fe q)m—a-l (mt) (r) dr, (3.2) 
= -L 

It has been demonstrated that this modified fractional-order derivative possesses two im- 
portant properties: the preservation of periodicity and the short memory, which considerably 
reduces the cost of numerical computations. Furthermore, it has been proven that contrarily 
to fractional autonomous systems defined using classical fractional derivative, the fractional 
autonomous systems in terms of the modified fractional derivative can generate exact periodic 
solutions. 
In order to generalize this work, one introduces in this section a similar modification to both 
Caputo fractional-order derivative and Riemann-Liouville fractional-order derivative as fol- 
lows. 


Definition 3.3. (The Caputo fractional derivative with sliding fixed memory length) Let a > 0, 
L > 0, m an integer such that m = [a] +1 and f € C™"{a—L,b]. The Caputo fractional 
derivative with sliding fixed memory length is defined by 


1 


(ma), 


aay ee ade. (3.3) 





MCDEF(H) = = 


bo 


Definition 3.4. (The Riemann-Liouville fractional derivative with sliding fixed memory length) 
Let a > 0, L > 0, m an integer such that m —1 < a < mand f is a continuous function in 
[a — L,b|,the Riemann-Liouville fractional derivative with sliding fixed memory length is de- 
fined by 
MRL, 1 a” i. m—a—1 
eee — d A 
LDS) = panaam [to Flee, (3.4) 
Remark 3.5. From (3.2) and (3.3) one gets 
m—1 OG _ Lye 


“De = Me De 
L f(t) L f(t) Xu T(k-—a+1) 
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Proposition 3.6. Under the assumption that the function f(t) is m-times continuously differentiable 


MRL 


LDri= 0 





SS Ail ae a 


=n T(k—a+1) ° 


Proof. By differentiation and performing repeatedly integration by parts, one has 





"L Dif (t) = Fm a) D a fe —1)"-"-1 ¢(7)dr, 
(m—1) pm—a-1 a iF qm-l t 
— Le "I ; , T(m— a —1) den—1 mc me (0 a 





m—1 k—a : 
re es + Ew [ @-2 far, 


setting I = Ter (tt) ~*!f(r)dt, and performing successive integrations by parts 


one obtains 








_f=LE* 1 t ae 
"= Ta-a) Ta =a [ @-0 fi'(t)dt, 


5 ft=L)L> FGshE 1 ff (tt) 1) (a) dz, 











T(1—«) rea) TO =o) Jee 
y yi = os : : —a+m-—1 ¢(m) 
ry. Tike) | T'(m — «) | @-0 ff" (a) dt, 
="" Di f(t), 


Therefore 





Remark 3.7. From (3.5) and (3.6) one has 


MRL 


Dis PHO = , DAO 





= f (f= L)LE-« (3.7) 


= teoes i) ~ 
In the following parts, one denotes the operators of Caputo and Riemann-Liouville frac- 
tional derivative with sliding fixed memory length by /D*. 


3.1 Fractional derivative of some elementary functions 


In order to highlight the amazing properties of the fractional derivative with sliding fixed 
memory length one consider two elementary functions (the power and exponential functions), 
for which one computes their new derivatives. 
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3.1.1 New fractional derivative of the power function 


Let f(t) = t", n € IN*,« >0,L > 0 and m is an integer such that m—1 <a < m. 
If n < m, then f(")(t) = 0, substituting in (3.3) yields MD*(t") = 0. 
If n > m then by repeated integrations by parts of the relation (3.3) one obtains 


n—-m aN am aia = Lys 


Mpa(yny) — 
a age Dae m—k)IT(—a+m+k+1)’ 





(3.8) 





Remark 3.8. (Fractional derivative of a constant function) 
If f is a constant function (i.e. f(t) = C for all t € [a — L,b], and C any constant including zero) 
then one has 

VO =e. 


3.1.2 Fractional derivative of the exponential function 


Let f(t) =e = » f «>0,L>Oand mis an integer such that m—1 <a <m. 
pl & 
70 








co p—m Loatmt+k(- — [)\pom—k 
MDK(e) = wt) 
pad ko (P m—k)IT(-«+m+1+k) 
p-m 


foe) Lees —s Lier’ 

7 ud (p—m—k)'T(k-a+m+1)’ 
3 y Leer e -_ L)P-* 

((p=k) k= ae)" 


co pa +m (t — L)P | 0° magmas ees 2) 
7 -—. ga PIC (—# + m +2) a 


-(¥ —L)? 2 Loatm+k 
pao PIT ( rainy A Sense Tee ey 
Lk 
eb a+m 
ar (—w+m+1+k)’ 
= ELE a(L) 

















a 





3.2 Derivative of a periodic function 
The main result of this paper is stated in the following theorem. 


Theorem 3.9. Let « > 0, L > 0 and m an integer such thatm—1<a«< mand f € C™[a—L,b]. 
If f is a periodic function with period T, Then MD¢ f is a periodic function with the same period T. 


Proof. Suppose that f is a periodic function with a period T. The aim of this proof is to 
demonstrate that the function g(t) = “D*f is a periodic function with the same period T (ie. 
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g(t+T) =g(t)). 
One has 


t+T 
s+) = MDA f+ TD = remy fT em et Tae 


= Fm nw) = [ (b= sym 1 f0 (8 + 27s, 
[91 f0 (as, 


= MDff(t) = g(t). 


Thus, “D*f is a periodic function with the same period T. 


3.2.1 Fractional derivative of some fundamental periodic functions 


Note first that the functions “°D* sin(t) and “Df cos(t) have been calculated in [1]. 


Example 3.10. (Fractional derivative with sliding fixed memory length of the sine function) 
By definition 





oi FYE _ LL 


MDE f(t) = “LDEF(E) Tkoa bi) 


k=0 


Therefore 


mal a ( (sin(t — L))L*-# 
Mp* sin(t) = “"; Df sin(t) i ag )) 























— —a+1) ’ 
= L~“ sin(t — Bia ne oe eee 
[24] (—L2)k [*] (—L?) 
L-* imij=Lj—1s* t—L 
dX FOk+1—a) sin( ) u T(@k+2—2) cos( ), 
- » (El (oe 
= L “sin(t — L)(E21~-4(—L*) Xu rorai— a)? 
ie 5 ["3] (—L?)* 
+L’ cos(t — L)(E22-«(—L*) Xu TQk+2—a) ) 
= asin(t — L) + bcos(t — L), (3.9) 
wet] (_2)k ["r]  (—12)k 








where, a = L~*(E1-4(—L?) — © \e a TP Fog pO) 


», T@k+1—a) &, T@k+2—a)” 
One notices that, “'D# sin(t) is a periodic function with the period T = 27t. This analytical 
result is displayed in figure (3.1), for some values of a and L = 32.1. 


Example 3.11. (Fractional derivative of cosine function) 
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Figure 3.1: Fractional derivative of the Sine function for L = 32.1 and some values of a. 


By definition 


1 e (tf — L))Lk-# 
m= dik COs 


=. Teorey l) 
= L* cos(t — L)Eo1_,4(—L”) — L1* sin(t — L)E22-«(—L’) 


MDE cost) = wD cos(t) 






































ee, ale) (aye 
af X foi u TQk+2—a) nt 2), 
. ») 2! ny 
=L * cos(t= L)(Eat-2(—L") Xu Tee l= ae 
1-a g: 2 ea (Gi 
— L* “sin(t — L)(Ex2-4(—L*) Xu T(2k+2—a) ) 
= acos(t— L) —bsin(t — L), (3.10) 
where, 
[= 2)k 
a =L™"(Ep;-2(—-L*) Xu ae — x) y 
and 
eee ~ (=F , 
& Tk +2—a) 


Obviously “Df? cos(t) is a periodic function with period T = 27. 


3.3 An interpolation property 


It is known that the operator of Griinwald-Letnikov fractional derivative with sliding fixed 
memory length is an extension of the integer-order operator 47, (see [1]). 
The following proposition proves that the Caputo and Riemann-Liouville operators of the 
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fractional derivative with sliding fixed memory length verifies this property for « — m, but 
not for a — m—1. 


Proposition 3.12. Let L > 0 and0 <m-—1<.« < msuch that m is an integer number, and let f(t) 
having (m +1) continuous bounded derivatives in [a — L,b]. Then , for all t € [a,b], one has 


lim MDEF(t) = f(t), 


a—->m 


and 


Proof. One has 


lim MD¢F(t) = lim Fa) rh (t—r)"-*1¢™ (x), 


am asm T t-L 


Ry re ies Ae), 
Se Te aD) 


f t-or ep (aan, 


t—L 
=f E-L+ fpr @ar, 
= 


+ i 
am Tm — a +1) 





For a + m—1, one has 


li Me ee (t) = li a | G22 (") (t)dt 
Pastas L if pn ee) t f 


Example 3.13. 
Let f(t) =e’, then 
Dee ee ea og a 


Therefore, 
lim VD Se (es =p. 
However, 


im, Mpte Se “LE o(L) =e Fel = 1), 


=e! — ef b= fl) (ty) — fDi — 1). 


Example 3.14. 
Let f(t) = #", one has 


n—m Ol mapas = Hew 


"D} (t") = u (n—m —k)T(—a+m+k+1)° 
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Putting N =n—mandt—L =a, then 


a 7 N niLka N-k 
jim 7D} (t") = ay (N— bik’ 
n! SN NILkaN-k 
N! a (N — k)tk!’ 























_ N n! wai 
= 7 (4 + £) = ar 
m 
= Fras) 
dt 
However, 
N {pk+1gN—-k 
Fi Mpnasyn\ nN. a 
ee )= d= BNREDY 
nl x (N +1)!LkaN+1-k n! (t—Ly 
(N41)! & (N+1-4)Ik! (n—m-+1)! : 
es nt ~(m-1) 
“Weil  G@-=m—-pi- 
_ n! n—(m—1) ni n—(m—1) 
“Gam a= 
qm-1 qm-1 





= tt = (tb = fH) = f= L). 


3.4 Comparison between some results of classical fractional-order derivatives and 
fractional order derivatives with sliding fixed memory length 


The previous results are summarized in the table (3.1), in order to highlight the differences 
between classical fractional-order derivative and fractional-order derivative with sliding fixed 
memory length. 


3.5 Fractional-order autonomous system with exact periodic solution 


As previously mentioned, any autonomous fractional-order system expressed in terms of clas- 
sical fractional derivatives cannot have any exact periodic solutions [13, 22,23]. 

Conversely to these results, one presents some examples showing that fractional-order au- 
tonomous systems (linear and nonlinear) expressed in terms of fractional derivatives with 
sliding fixed memory length can have exact periodic solutions. 


Example 3.15. (Linear fractional-order system) 
Let consider the following linear fractional-order autonomous system 


MDX (t) = AX(t), (3.11) 


a 


b 





= [x1] 
where X(t) € IR and A = ( i ) , with a = L~*(E)_-4(—L*) y T Cy ), 
p=0 





m—2 ] 


ay Z —L?)P 
b= LI *(Ezaa(-E?)— ES rpaaa) 
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Classical fractional derivative Fractional derivative with sliding fixed 
C D* or aoe memory length “Di 
(k) ( « 
i= De pi eu s “Dtf() = “DEF 
fim RUDFF(®) = lim SDsf( = FW) lim PDFF® = FO) 
: RL pa m—1) : Mpa 
im a D F(t )= f! (t), am, 7 Di f(t) 
lim gDEf(t) =f") (E) — fF" @) SPN) ape Mest) 
a->m— 
T 1 ae n—m {L-atm+k mG n—m—k 
NDE") = GDP(E") = rye : ee) u, = aeraererasy 
oe = = Tea (t = a) # 0, 
CDi =U Mie =o 
ALD? sint = f' *Eps-y(—t*) MD* sin t = asin(t — L) + bcos(t — L). 














Table 3.1: Comparison between some results of classical fractional-order derivatives and frac- 
tional order derivatives with sliding fixed memory length. 


t 
- For L = 2k7, where k is a non-zero integer. The vector function X(t) =c ( me ) pcER 


is an exact 271—periodic solution for the system (3.11). 
By definition, 


we eostt) 
M a 
DIXG) Se) ; 
MDExt) =e ( Hp) ) 


Then, from (3.9) and (3.10) one obtains 


ee 7 acos(t — 2k7t) — bsin(t — 2kzr) 
eC ( asin(t — 2k7t) + bcos(t — 2k7t) ]’ 


(he) (Sat ate): 
=A (Sanat) 
= AX(t). 


Therefore, X(t) =c ( coos) 
sin(t) 


- For L = 5, one has 


) is an exact 27t—periodic solution of (3.11) with L = 2kzc. 


MDex = pa eae a Hay 
2 


asin(t — ) + bcos(t — 


( 

a asin(t Best ), 
( -. 
B 


—acos(t) + bsin(t) 
(<4) Cane ): 


=6( on 
= BX(t), 
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‘ - b a cos(t) 
with B = ( ae a sin(?} 
exact 27t—periodic solution of the system 7 DEX(t) =DX(e), 


) # A. Thus, X(t) = ( ) is not solution of (3.11), but it is an 


Example 3.16. (The predator-prey model with Holling type II response function) 

All population species possess the property of heredity, which means the passing on traits 
from parents to their offspring, either through asexual reproduction or sexual reproduction. 
The offspring cells or organisms acquire the genetic information of their parents through 
heredity. This property makes fractional differential systems models more efficiently regard- 
ing some specific problems than ordinary differential ones. 

Motivated by this fact, we introduce the fractional version of the Holling-Tanner model [21] 





as follows 
xy 
el 8) ate (3.12) 
D*y = ray(1— <2). 


Where D: denotes a standard fractional-order derivative operator and « € [0,1] is the fractional- 
order related to the hereditary property of the population (a value of w close to an integer 
number means that the population has a weak hereditary property), x(t) > 0 and y(t) > Oare 
the density of prey and predator populations at time ft respectively. The parameters 1; and rz 
are the intrinsic growth rates, K represents the carrying capacity of the prey, q is the maximum 
number of preys that can be eaten per predator per unit of time, m is the saturation value (it 
corresponds to the number of preys necessary to achieve one half the maximum rate q), y is a 
measure of the quality of the prey as a portion of food for the predator. 


Since exact analytical resolution of this nonlinear system is unavailable, one resorts to 
qualitative and numerical study. For this purpose the parameters are set to 7) = 1, 12 = 
02; K- = 25, 9 = g, m = 1and y = 0.95, the system (3.12) has two equilibrium points 
Eo = (25,0) and E; ~ (7.1429, 6.7857). 


¢ The characteristic polynomial of the Jacobian matrix evaluated at Eo is given by 
P(A) = A* +ayA + a2 = A* +0.8A — 0.2. 
So a2 = —0.2 < 0, then according to Proposition 1 in [7] Ep is unstable for all « € [0,2). 
e The characteristic polynomial of the Jacobian matrix evaluated at E; is given by 
P(A) = A* — 0.1409A + 0.0747. 


So a, © —0.1409 and a2 & 0.0747 > 0. 

Applying Hopf-Like Bifurcation theory [2-4] and using Proposition 1 in [7], one obtains 
the Hopf-Like bifurcation value 

2 /a2 
The fixed point Ej losses its stability, and a periodic motion (S—asymptotically periodic 


for the classical fractional derivative and exact periodic for fractional derivative with 
sliding fixed memory length) appears. 


2 
a* = —cos*( 
7 





) © 0.8341. 
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(c) (d) 
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AL J 
§ [ae nee eS crer 

aL J 
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(e) (f) 


Figure 3.2: Time evolution and phase portrait of system (3.12) for a = 0.9 (a,b) 
S—asymptotically T—periodic solution with T * 27.2 for classical fractional operator. (c,d) 
Exact T—periodic solution for the fractional derivative operator with sliding fixed memory 
length. (e,f) Comparison between the two solutions. 


To illustrate these results, one solves the system (3.12) numerically by developing a Matlab 
code using a discretization technique based on the formula (3.7). 

Choosing a value for a greater than a*, for example, « = 0.9, one compares the solution 
of (3.12) in terms of classical fractional operator and its solution in terms of the fractional 
operator with sliding fixed memory length L = 30. The two trajectories are start from the 
same initial point Xp = (2.64,4.88), belonging to the attracting limit cycle. The results are 
shown in Fig. 3.2. 

An S—asymptotically T—periodic solution with T ~ 27.2 is obtained for classical fractional 
operator as shown in Fig. 3.2(a,b), and an exact T—periodic solution is obtained for the 
fractional derivative operator with sliding fixed memory length as shown in Fig. 3.2(c,d). 


4 Conclusion 


In this article, one modifies the Caputo and Rieman-Liouville fractional-order derivatives by 
fixing the memory length and varying the lower terminal of the derivative. It is shown that 
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the modified fractional derivative operator preserves the periodicity. Consequently, periodic 
solutions can be obtained in fractional-order systems expressed in terms of the new operator. 
Two examples are investigated to highlight this property for a linear system provides an 
analytic expression of an exact periodic solution is computed and for another nonlinear system 
for which exact periodic solutions are obtained using qualitative and numerical methods. 
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